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To date, quantum computational algorithms have operated on a superposition of 
all basis states of a quantum system. Typically, this is because it is assumed that 
some function /is known and implementable as a unitary evolution. However, 
what if only some points of the function /are known? It then becomes important to 
be able to encode only the knowledge that we have about/. This paper presents an 
algorithm that requires a polynomial number of elementary operations for 
initializing a quantum system to represent only the m known points of a function/. 

PACS numbers: 03.67 .Lx, 03.67.-a, 89.80.+h, 03.65.Bz 



The study of computation as a physical process has produced fascinating results in the form of 
new approaches to problem solving and information processing that exhibit impressive speed-up 
over classical approaches for some problems [1—6]. All of these algorithms operate on a 
superposition of all basis states, usually assuming prior knowledge of some function/that can be 
effected as a unitary evolution of the system. However, it is often the case in computational 
problem solving that a function /is not known but rather that some example points of the function 
are all the knowledge available. Such situations are particularly common in the field of 
computational learning, which typically deals with learning a function from examples. If quantum 
computation is to become a broadly applicable paradigm, it must be possible for quantum 
algorithms to function given only a limited set of functional points. The first major step toward 
this goal is developing the ability to encode a set of examples of a function as the state of a 
quantum system. If this can be achieved, then the possibility of quantum computational learning 
algorithms presents itself. For example, fourier-based machine learning algorithms such as those 
by Kushilevitz and Mansour [7] and Jackson [8] are tempting candidates for generalization into the 
quantum realm. In fact, Bshouty and Jackson investigated a quantum computational approach to 
learning [9]. However, their algorithm lacks the capability proposed here (relying instead upon a 
quantum example oracle that somehow encodes the entire function/even though the function is not 
known a priori) and therefore, though theoretically interesting, is not a plausible quantum 
algorithm. This paper presents an algorithm for encoding, in polynomial time, a set of examples of 
a function as the state of a quantum system. 

Given a set Tof m examples of a function/ the goal is to produce 

/)= Lf(z)\z) (1) 

as the quantum state of n qubits. It will be shown that the state f^j can be constructed using a 
polynomial number (in n and m) of elementary operations on one, two, or three qubits. In what 
follows, the specific qubits to which an operator is to be applied are indicated as subscripts on that 
operator. For simplicity, consideration will first be restricted to the case of /:z — >s with 
z g {0,1 }" and ss {-1,1}. First, define the set of 2-qubit operators, 
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where se {-1,1 } and m<p <1. These operators form a set of conditional Hadamard-like transforms 
that will be used to incorporate the example set into a coherent quantum state. There will be a 
different S s ' p operator associated with each example in the set T. Next, define 

"0 1" 

1 0_ 

which flips the state of a qubit. Now, define the 2-qubit Control-NOT operator 

F 
/ 2 

where and I 2 

second qubit if the first qubit is in the |0) state. Similarly define F l which conditionally flips the 
second qubit if the first qubit is in the |1) state [the matrix representing F l is the same as Eq. (4) 
with I 2 and F exchanged]. Finally, introduce four 3-qubit operators (really just four versions of 
the Fredkin gate [10]). These operators are used to identify specific states in a superposition, 
similar to Grover's identification of which state should be phase-rotated in his search algorithm 
[2] . The first of these, 
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t 2 are the 2x2 zero and identity matrices respectively, which conditionally flips the 
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where the 6 are 6x2 and 2x6 zero matrices and I 6 is the 6x6 identity matrix, conditionally flips 
the third bit if and only if the first two are in the state |00). Note that this is actually equivalent to 

^ 00 ^ 00 

F , where F is the 3-qubit generalization of Eq. (4); however, this operator will be used with 
the third qubit always in the |0) state and thus can be thought of as performing a logical AND of the 
negation of the first two qubits, setting the third to |1) if and only if the first two are |00). The 

~ A | A 1A ^11 ^00 ^ 

other 3-qubit operators, A , A , and A , are variations of A in which F occurs in the other 
three possible locations along the main diagonal. 

The algorithm for encoding Tinto a quantum system requires n + (n - 1) + 2 qubits, arranged 
in three quantum registers labeled x, g, and c. The quantum state of all three registers together is 
represented as \x,g,c), and the algorithm proceeds as in Fig. 1. 
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FIG. 1. Initializing amplitude 
distribution in a quantum system 



The x register will hold a superposition of the examples in the set T— there are n qubits in the 
register, and the value of / for an example will be used as the coefficient for the state 
corresponding to that example. The g and c registers contain only ancillary qubits and are restored 
to the state |0^ by the end of the algorithm. A high-level intuitive description of the algorithm is as 
follows. The system is initially in the state |0^. The qubits in the x register are conditionally 
flipped so that their states correspond to the first example. The S s ' p operator corresponding to that 
example then changes the basis of the system such that one of the states has a coefficient that 
matches the example's value for/. This same state then is changed to one that exists outside of the 
subspace affected by the S s ' p operators, in effect making it permanent, and the process is repeated 
for each example. When all the examples have been processed, the result is a coherent 
superposition of states corresponding to the sample points, where the amplitudes of the states all 
have the same magnitude but have different phases according to the values of their corresponding 
examples. This is somewhat reminiscent of the first part of Bernstein and Vazirani's fourier 
sampling mechanism [11]. However, they assume that the function/is known and calculable in 
polynomial time, whereas here we do not know the function/, but only a small set of examples 
drawn from/ 

In analyzing the complexity of the algorithm it is assumed that line 1 can be done trivially. The 
loop of line 2 is repeated m times and consists of n (lines 3-5) + 3 (6-8) + n-2 (9-10) + 1 (11) + n- 
2 (12-13) + 1 (14) operations, and line 15 requires one more operation. Thus, the entire algorithm 
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requires m(n + 3 + n-2 + 1 + n-2 + 1) + 1 = m(3n+l)+l operations and is O(mn). This is optimal 
in the sense that just reading each instance once cannot be done in any fewer than mn steps. 

A concrete example of a simple 2-input function will help clarify the preceding discussion. 
Suppose that we are given the set T= {/(01)=-1,/(10)=1,/(11)=-1}. The initial state is |0^, and 
the algorithm evolves the quantum state through the series of unitary operations described in Fig. 
1. For convenience lines 3-6, 8-10, and 12-14 of the algorithm are agglomerated as the compound 
operators FLIP, AND, andANDt respectively and are treated in more detail in Eq. (6-8) and Figs. 
2 and 3. First the qubit states of the x register are flipped to match the first example, and the c 
register is marked so the state will be affected by the S s,p operator. 

1 00, 0,00) FLIP > 1 01, 0,10) 
Next, the appropriate S s ' p operator (with s equal to the output class of the instance being 
processed, here -1, and p equal to the number of the instances, including the current one, yet to be 
processed, here 3) is applied, creating a new state in the superposition. This process will be 
referred to as state generation and corresponds to line 7 of Fig. 2. 

S — ) -^|01,0,ll) + |01,0,10) 
The two states just affected by the S s ' p operator are marked in their g registers. 

AND > -^|0Ull) + j||0U10) 
One of the marked states is made permanent by setting its c register to 1 01) and this state now 
represents the first example. The c register of the other state is returned to |00), and it is ready to 
generate a new state. This is performed by line 1 1 of the algorithm. 

F ^ ) -^|01,l,01) + ||01,l,00) 

Finally the work done in the g register is undone and at this point one pass through the loop of line 
2 of the algorithm has been performed. 

ANDt ) -^|01,0,01)+ J01,0,00) 

Now, the entire process is repeated for the second example. Again, the x register of the 
appropriate state (that state whose C2 qubit is in the |0) state) is selectively flipped to match the new 
example, and its c register is again marked. This time the selective qubit state flipping occurs for 
those qubits that correspond to bits in which the first and second examples differ — both in this 
case. Notice that the state representing the first example is not affected by this operation, nor any 
of those that follow. 

FLIF > -^101,0,01) + 1|10,0,10) 
Next, another S s,p operator is applied, this time with s = 1 and p = 2. 

— 2% -^101,0,01) + ^ j||io,o,ii)+jrjf|io, o,io) 

The two states just affected by the S L operator are marked in their g registers. 

AND ) -^|01,0,01) + ^|10,l,ll) + Ji|10,l,10) 
The c register is used to save another state (this one representing the second example). 

^ ) -^|01,0,01) + ^|10,l,01) + ||10,l,00) 
A final bit of house keeping resets the g register, preparing for the process to be repeated again. 
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^|01,0,01> + ^|10,0,01) + #0,0,00) 



Finally, the third example is considered. The x register of the generator state is again selectively 
flipped and its c register is marked. This time, only those qubits corresponding to bits that differ 
in the second and third examples are flipped, in this case just qubit X2- 

FUP > -^|01,0,01) + ^|10,0,01) + ||ll,0,10) 

Again a new state is generated to represent this third example. Notice that the generator state is 
now left with a magnitude of 0, because there are no more examples to process. 



^|01,0,01) + ^|10,0,01>-^^|ll,0,ll) + J°jl|ll,0,10) 



Once again the house keeping steps are performed 
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Finally, line 15 performs one last operation, restoring all the ancillary qubits to their initial state. 

Fc2 ) -^[01,0,00) + ^[10,0,00) -^|1 1,0,00) 
At this point, the g and c registers are not entangled with the x register, and therefore the system 
can be simplified as 

-^|01) + ^|10)-^|11), 
and it may be seen that the partial function defined by the set Tis now represented as a quantum 
superposition in the x register. 

Here the complex operators FLIP, AND, and AND^ are specified in more detail. In operator 
notation, 
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(6) 



Alternatively, FLIP may be represented as a quantum network as in [12]. Figure 2 shows one 
such network for a 3-input function. In the figure, a '?' indicates that the F operator is applied 
only if the value of the bit in question differs from the value of that bit for the previous example. 
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FIG. 2. Quantum network for 
implementing FLIP 



Also, in operator notation 
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AND and AND f also may be represented as quantum networks, and corresponding example 
diagrams appear in Figs. 3(a) and 3(b) respectively. Note that the only difference between the two 
is the order in which the qubits are operated upon since all the A operators are their own adjoints. 




FIG. 3. Quantum networks for 
implementing AND and AND^ 

The quantum network representation of the entire algorithm requires the addition of gates for 
performing the S and F operations, and the complete network is shown in Fig. 4. The network 
shown will be repeated several times — once for each instance in the training set. Each such 
repetition will make use of different A operators and flip different qubits using the F gates, 
according to the particular example being processed. 




FIG. 4. Quantum network for encoding a single 
function example 

The trick that makes the algorithm work is the fact that whenever an S s ' p operator is applied, 
there are no states with the c register in the |1 1} state. This is crucial to the construction of f^j 
because it allows the operations to be unitary without generating "extraneous" states. Another key 
is the ability to identify a specific state in the superposition (the one just generated) using the AND 
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and AND f operators; this provides the ability to mark that particular state and only that state so as 
not to operate on it again. And this is done reversibly so that no disruptive entanglement results. 

The algorithm of Fig. 1 can handle only binary functions. However, generalization to 
functions with N input/output values is straightforward. The only necessary changes are a 
generalization of the x register (and only the x register) to include systems with more than two 
states, generalization of some of the operators and the representation of more than two output 
values as points on the unit circle in the complex plane (representing a binary output with the 
values -1 and 1 is just a special case of this). The operator generalizations are as follows. 
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where now s& {0, AM } and m<p <1. The unit circle in the complex plain is partitioned into N 
pie pieces, and the points on the pie piece borders correspond to the function values. By 
convention, one of the values is always on the positive real line. 
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(10) 



where 6 and I N are the NxN zero and identity matrices respectively, and C is the NxN circulant 
matrix obtained by shifting all rows of I N down one, with the Mh row becoming the first. There 
is now a set of F operators, one for 2-state qubits, one for 3-state qubits, etc., and the F 
operator defined in Eq. (4) is the special case for N=2. These operators act on a pair of qubits, the 
first of which is in the c register and thus is still just a 2-state system, and the second of which is in 
the x register and thus now has N states, as N F°\0q) = \0(q + l)modN) and N F°\lq) = \lq), 
conditionally rotating the second qubit (labeled q) to the next highest value (but rotating the 
maximum value back to 0) if the first qubit is in the |0) state. Note that F l does not need to be 
generalized since it only acts on qubits in the g and c registers, which are still only 2-state qubits. 
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where a = 2N 2 -2, the 6 are 2xa and ax2 zero matrices, I a is the axa identity matrix, F is as 

2 2 ■ ^ 00 

defined in Eq. (3) and the entire matrix is 2N x2N . Again there is now a set of A operators 
and the A operator defined in Eq. (5) is the special case for N=2. These operators act on three 
qubits, as N A°°\00r) = \00(r + l)mod2> and N A 00 \ pqr) = \ pqr) for p,q * 0, conditionally 
flipping the third qubit (labeled r), which will still be a 2-state qubit as it is in the g register, if the 
first two qubits are in the 1 00) state. Note that whereas the operator A 00 defined in Eq. (5) was 
really equivalent to a 3-qubit version of Eq. (4), this is not the case with these generalized 
operators. In other words, Eq. (1 1) is not equivalent to a 3-qubit generalization of Eq. (10). This 
is because N F° now affects qubits (in the x register) with more than 2 states, therefore requiring a 
rotation through the various states rather than a simple flip between 2 possible states; whereas 
N A 00 still affects qubits with only 2 states and thus simply toggles between states. Since this 
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operator will be used with the third bit always in the |0) state, it can still be thought of as 
performing a sort of logical AND (actually a generalized AND where the inputs can be greater than 
binary) on the values of the first two qubits, setting the third to |1) if and only if the first two are 
1 00). The other A operators must also be generalized, a set for each possible pair of values that 
any of the qubits in the x and g registers can assume. In other words we need all the sets of 
operators N A^ for 0<iJ<N-l. However, these sets of operators are simple variations of the N A°° 
set given in Eq. (11) with the sub-matrix F occurring in different locations along the diagonal 
according to the values for i and j . 

The algorithm for the more general case is the same as in Fig. 1, except for two changes. 
First, the operators must be replaced with their generalizations defined in Eq. (9-11), and second, 
the IF statement of line 4 must be changed to a WHILE statement. This is again because instead of 
simply flipping a qubit's state between two possibilities as the binary operator F defined in Eq. 

at n 

(4) does, the generalized operator F defined in Eq. (10) rotates a qubit's state through N 
possibilities in a specific order one at a time. Thus, the compound FLIP operation defined by lines 
3-6 of the algorithm now may require up to (N-l)n operations instead of the 0(ri) required for the 
simpler binary case. This brings the total number of operations required for the algorithm up to 
(N+l)mn + m + 1, which is still OimnN). This reduces to Oimn) assuming N«m,n . 

In summary, this paper presents a polynomial-time algorithm for initializing a quantum system 
to represent m known points of a function/. The result is a quantum superposition with m non- 
zero coefficients — the creation of which is a nontrivial task compared to creating a superposition of 
all basis states. It may be appropriate to mention here that very recently work as been done to 
analyze Grover's algorithm [2] for the case of arbitrary initial amplitude distributions [13]. 
Further, the paper suggests a new field to which quantum computation may be applied to 
advantage — that of computational learning. In fact, it is the authors' opinion that this application 
of quantum computation will, in general, demonstrate much greater returns than its application to 
more traditional computational tasks (though Shor's algorithm is an obvious exception). We make 
this conjecture because results in both quantum computation and computational learning are by 
nature probabilistic and inexact, whereas most traditional computational tasks require precise and 
deterministic outcomes. 
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